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ABSTRACT 

Extinction due to intergalactic grey dust has been proposed as an alternative to accelerated expansion to account for the dimming of 
Snla fluxes beyond z — 0.5. The "replenishing" grey dust model, although fitting the observational data, does not seem to be based on 
physical assumptions. For this reason, in this paper we propose a new grey dust model whose dust distribution follows the comoving 
SFR density evolution, a reasonably established phenomenon. This new model is compared with the updated photometric data sample 
from the High Z Supernova Search Team (HZT). Also, using pairs of supernovae at different redshifts, the possibility of any "patchy" 
distribution proposed to explain the discrepancies between the "high z" dust model and the Snla observations at very high redshift 
(z > 0.9) is ruled out. Finally, the data are compared with different models with and without dark energy and a best fit to a universe 
with adimensional parameters fi m0 = 0.31 and Q. A0 = 0.69 is obtained, from which an age of the universe of t 14.6 x 10 9 years 
is derived. This age is compatible with the age of the globular clusters using an equation of state oj = - 1 and it obviates the need to 
resort to any kind of phantom energy. 
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1. Introduction 

After the discovery by Riess et al.(1998) and Perlmutter et 
al.(1999) of the excessive weakness of type la supernovae (Snla 
) at high redshift, cosmology has undergone a revolution simi- 
lar to the cosmic microwave background (CMB) discovery. The 
cosmological interpretation of these observations led to the pos- 
tulation of the accelerated expansion of the universe as a con- 
sequence of a "dark energy" of negative pressure and unknown 
origin. Given the difficult physical interpretation of this dark en- 
ergy, either in terms of a cosmological constant, quintessence, 
or even some form of phantom energy, several mechanisms have 
been proposed to explain Snla observations without resorting to 
an accelerated expansion. Among them are models of intergalac- 
tic grey dust, which try to explain the dimming of the observed 
brightness at z — 0.5 by assuming that an absorbent medium ex- 
ists in intergalactic space whose extinction does not depend on 
wavelength. These dust models have been driven to evolve at the 
same rate as the data supplied by Snla teams. 

The first models, proposed by Aguirre (1999a, 1999b), were 
based on both a model of dust formed by carbon needles and 
a Draine & Lee dust model (Draine & Lee 1984), obtaining a 
density parameter of Qoust ~ 10~ 5 . Also, Goobar, Bergstrom 
& Mortsell (2002) have developed a model of intergalactic dust 
whose effects on the distance modulus vs. z are indistinguishable 
from those caused by a ACDM universe with £2 m o — 0.3, Q.\o =s 
0.7 (more information on the intergalactic medium amount and 
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its effects can be found in, for example, Bianchi & Ferrara 2005; 
Inoue & Kamaya 2003; Nath, Sethi & Shchekinov 1999). 

In any case, grey dust models imply universes with Q m o = 1 
or non-Euclidean geometries. The first possibility can be dis- 
carded by the 2dF experiment (Peacock et al. 2001), while the 
second can be rejected by analysis of the anisotropies observed 
in the CMB (Spergel et al. 2003), from which a total parameter 
Q. = 1 is derived. 

Apart from grey dust models, another effect that could 
hamper the need for an accelerated expansion universe to ex- 
plain Snla data is their luminosity evolution. Observational evi- 
dence shows Snla luminosity variations with galactic morphol- 
ogy (Hamuy et al. 1996, 2000; Branch et al. 1996; Reindl et al. 
2005), while some theoretical studies (Hoflich et al. 1998, 2000) 
suggest that Snla luminosity could have some dependence on 
redshift. Were this true, this would question the status of stan- 
dard candle for this kind of objects. In any case, these studies 
are beyond the scope of this article and are merely mentioned as 
further possibilities in addition to grey dust models in order to 
avoid accelerated expansion. 

In this paper we present observational evidence to refute the 
high-z grey dust model even in the case of patchy distribution. 
We also propose a new grey dust model with a density distribu- 
tion following that of the star formation rate (SFR) density. This 
SFR-based distribution is based on observational data, thus mak- 
ing the grey dust model based on physical assumptions that best 
fits the supernova data. In Sect. 2 we give an overview of ex- 
isting theoretical grey dust models, including high-z and the re- 
plenishing grey dust. In this section the new SFR-based distribu- 
tion model is also developed . In Sect. 3 we derive the dust-free 
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model by best-fitting the Snla data, obtaining its density param- 
eters and the age of the universe, and compare the data with dust 
models, including the one developed here. We also explore the 
possibility of reconciling the observational data with the high-z 
dust model by means of a patchy dust distribution. Finally, we 
draw our conclusions in Sect. 4. 



2. Theoretical models 

The existence of grey dust would have an enormous impact on 
the analysis of the High Z Supernova Search Team (HZT) and 
the Supernova Cosmology Project (SCP) observations, return- 
ing the standard model of the universe to the point where it was 
before 1998 with their contradictions of the CMB anisotropy ex- 
periments, which point to a flat universe (Q = 1). This is the 
reason a detailed study of the different possibilities is required. 

Any intergalactic dust model producing grey extinction rel- 
evant to the study of Snla observations must have the following 
characteristics: 

1 . Large grain size, to account for the achromaticy of the ab- 
sortion and to comply with the mechanism of preferential 
destruction of smaller size grains (Aguirre 1999b). 

2. An efficient expulsion mechanism in the galaxy where the 
grains are formed. It must explain why only this type of dust 
is found in the IGM. 

3. Coherence with the expected metallicity of the intergalactic 
medium. 

4. The dust distribution must be based on a hypothesis in turn 
based on reasonably established physical observables. 

2.1. High-z grey dust 

Even though Aguirre's (1999a, 1999b) high-z dust model im- 
plies a uniform dust distribution similar to that of non-relativistic 
matter, 

Pdust=pL t (l+^ (1) 

which fulfils the first three of the four requirements mentioned 
above, the latest results published by the HZT (Riess et al. 2004), 
with Snla data beyond z = 1.5, almost certainly (at the Act level) 
reject the possibility of this high-z dust model. Figure 1 shows 
the large deviation of the Aguirre (1999a, 1999b) model in the 
high-z region compared to the experimental data and the best fit. 

However a patchy dust distribution could render this con- 
clusion unsound. If the intergalactic dust distribution were not 
homogeneous, it might be possible that the light coming from an 
Snla at high redshift (z > 0.9) propagates through more or less 
absorbent medium, depending on the source position in the sky. 
The deviation of the observed data from the high-z dust model 
would then be due to a bias produced by the grey dust, leav- 
ing the magnitude of the absorbed supernovae below the detec- 
tion limit of current telescopes. In this case, the sample observed 
at very high redshift would belong to the group of Snla whose 
optical path has not been very affected by the aforementioned 
absortion, and for this reason its apparent brightness would be 
higher than that inferred by supposing a completely homoge- 
neous dust distribution. 

The inhomogeneity in the grey dust distribution is a solid 
enough possibility to take it into account in spite of Aguirre's 
(1999a, 1999b) arguments. In his theoretical study, possible 



mechanisms for preferential ejection of grey dust from the galax- 
ies and the temporal scale on which it could have been homo- 
geneously distributed, are explained. They are partly corrobo- 
rated by the work of Renzini (1997) on the metallicity of the 
intracumular gas. Nevertheless, it is difficult to explain how a 
dust distribution formed by baryonic matter can escape from the 
baryonic large-scale structures observed in the universe. Aguirre 
himself (1999a) suggests that the study of systematic differences 
between the flux received from Snla hosted by field and cluster 
galaxies could distinguish between a universe in accelerated ex- 
pansion and an effect due to grey dust, respectively, by restricting 
implicitly the homogenity to the intracumular space. 

2.2. Replenishing grey dust 

As different surveys reached more and more distant redshifts, it 
became evident that the high-z dust model was unable to explain 
the data. However, the reasonable values obtained by Aguirre 
(1999a, 1999b) for the time scale of the expulsion mechanisms 
of his grain models have promoted the birth of new models that 
fit the data much better. A good example is the replenishing grey 
dust model, proposed by Goobar, Bergstrom & Mortsell (2002), 
whose match to the data is at least as good as the Q A + models. 

This replenishing grey dust model is based on a grey dust 
extinction that follows a density distribution of the form: 

Pdust=pLt( 1+ z) a . ( 2 ) 
where 

_/3 ifz<0.5 
fl_ \0 ifz>0.5. 

This model requires extremely efficient formation and expulsion 
of the dust in the range z > 0.5 so that the dilution due to the 
expansion of the universe is exactly compensated by the metal 
formation in galaxies and its further removal. However, this be- 
haviour is not consistent with the current knowledge of the over- 
all star formation episodes that took place in the universe. As a 
consequence, a physical explanation of this replenishing mech- 
anism is still lacking. Moreover, the hitherto unjustified sharp 
transition from a — to a — 3 represents a dramatic end of 
metal formation and expulsion from the galaxy without a transi- 
tion era. 

2.3. SFR-based model 

To avoid these features of the replenishing dust model, a grey 
dust model based on the star formation rate (SFR) density evo- 
lution is proposed here. It will be shown that the proposed model 
fits the existing Snla data equally well, while at the same time 
providing a dust production based on the current knowledge of 
star formation evolution in the universe. 

Following the results compiled by Taniguchi et al. (2005) on 
the comoving SFR density evolution (Figure 2), it can be ob- 
served that the comoving SFR density (p SFR ) is approximately 
constant in the range 1 < z < 2. In the region < z < 1, the SFR 
density may be reasonably represented by a simple relation with 
redshift, with lower SFR at recent times. For simpler description, 
the SFR density curve is then represented by these two regimes. 
For our purposes this simple parameterization is adequate, but 
see Oda & Totani 2005 for more complex parameterization. 

In this context, the comoving SFR density is represented by 

f 10 1+z if < z < 1 
Psfr a | cons tant if 1 < z < 2. 
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Fig. 1. Comoving p SFR as function of redshift. Figure taken from 
Taniguchi et al. (2005) 



3. Results 

3.1. Cosmological constraints 

In this work we have used the latest published HZT data (Riess 
et al. 2004; Strolger et al. 2004), since it represents the most 
significant sample with photometric values to date. This sample 
has greatly benefited from the results of the Great Observatories 
Origins Deep Survey (GOODS) project, carried out using the 
Advanced Camera for Surveys (ACS) on board the Hubble 
Space Telescope (HST). The Snla classification in the groups 
Gold, Silver, and Bronze used is explained in full detail by 
Strolger et al. (2004). The most certain classification is the Gold 
class, characterised by a clear spectroscopic confirmation of the 
supernova type. 

In this section, the observational data are fitted by a model 
composed of non-relativistic matter and dark energy, neglecting 
the contributions of radiation and relativistic matter, which do 
not play a significant role in the universe in the redshift range 
considered. The magnitude used for the fit is based on the dis- 
tance modulus, 



m-M = 51ogD L +25, 



(4) 



The evolution of the SFR proper density (p P SFR ) is 



pI fr (z) 



I 



a 3 -H 



1/2 



-dz, 



(3) 



where a is the scale factor at present and the value taken for Hq 
is 65 km s _1 Mpc~'. 

Nevertheless, although the dust density derives from the SFR 
evolution, not all the dust escapes from the galaxy to the IGM. 
At least a fraction of the metals produced will be retained in low 
mass stars that never reach the core collapse supernova phase 
and do not present strong stellar winds. An estimate of the re- 
tained metal fraction can be obtained from the Salpeter initial 
mass function (Salpeter 1955) from 0.1 M e to 8 M G . This SFR- 
based distribution will then include, as explained above, a transi- 
tion at z - 1 between the two regimes of maximum and decreas- 
ing SFR: 



where m and M are the apparent and absolute magnitudes, re- 
spectively, and D L is the luminosity distance in Mpc. Assuming 
a flat geometry, as derived from the spectra of the CMB 
anisotropies, we have obtained, via an x 2 criterion, the best fit 
with the parameters £2 m o = 0.31 and Q A o = 0.69, values simi- 
lar to those found in previous studies (Riess et al. 1998, 2004; 
Perlmutter et al. 1999), and consistent with those derived from 
2dF results (Peacock et al. 2001). 

These values are only indicative, given the fit dependence 
with the rejecting criterion. In Figure 2 we show the Hubble di- 
agram of the entire sample, together with the best fit found. 
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The constants p^ ustl , and p® ust2 , which include the estimation of 
metals retained in low mass stars, must be fitted in the region 
(0 < z < 1) to match the expected deviation at z = 0.5, while 
they are chosen to guarantee the continuity of the curve in the 
region z > 1 . 

Although the dust is produced by the star formation, the rel- 
evant physical magnitude to evaluate is the opacity, not the SFR. 
The resulting extinction is then delayed with respect to dust pro- 
duction owing to the time interval elapsed from production to 
insertion into the IGM. This delay can be estimated from as- 
sumptions for the velocity of the ejecta. 

In order to fulfil the above-mentioned requirements, this grey 
dust model must present grain characteristics similar to those 
proposed in previous models, reproducing the expected values 
for R v and absorptions (Goobar, Bergstrom & Mortsell 2002). 
Finally, the injection mechanism and the preferential destruction 
of small size grains will be assumed to be those discussed by 
Aguirre (1999b, 2001). 




Fig. 2. Hubble diagram. Distance moduli of observational Snla data vs 
z. Overplotted are the fi m0 = 0.31, Q A0 = 0.69 universe, which is the 
best fit for a flat cosmology obtained from^ 2 statistics. 



It is worth noting that these values of mass and dark-energy 
density parameters, assuming a Hubble constant of 65 km s _1 
Mpc -1 , yield an age of the universe around to - 14.6 x 10 9 years 
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as derived using Cepa (2004), with the constraint oj — - 1 for the 
equation of state of dark energy (i.e. a cosmological constant). 
This result is very nearly compatible with the globular cluster 
age of 12.9±2.9 Gy found by Carretta et al. (2000) and its time of 
formation, without resorting to a phantom energy (Cepa 2004), 
although the error margin is still too large to determine the equa- 
tion of state of the dark energy. 

3.2. The new dust model 

The Snla weakness at z - 0.5 and their excess brightness at 
z > 0.5, representing the transition from an accelerated to a de- 
celerated universe, are best appreciated by representing not the 
distance moduli but the distance moduli differences with respect 
to the empty Milne's universe in constant expansion. 

The grey dust models are represented as the distance moduli 
of a flat universe with £2 m o = 1, Qao = 0, with an atenuation 
obtained by 
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Am dust = -2.5 log \exp{(2 J Pdust(zMz)dz) j 



(5) 



where r(z) is the comoving distance travelled by the Snla light, 
Pdust(z) is given by equations (1), (2) or (3), and C is a constant 
used for dimensional purposes. 

In the case of the SFR-based dust model proposed in the 
present work (3), it will be necessary to include a delay between 
dust formation and its incorporation into the IGM to contribute 
to the observed extinction. To estimate the time interval in which 
the dust ejected could reach a homogeneous distribution in the 
intracluster space, we assume an expulsion velocity for the dust 
grains of around 1000 km s _1 , following Aguirre et al. (2001) 
and Shustov & Vibe (1995). A delay of Af » 3 x 10 9 years is 
obtained, implying that the homogeneity of the dust produced at 
z - 1 is reached around a redshift of 0.3 < z < 0.4. 

Figure 3 shows the different universe models, together with 
the weighted average data. From this figure it can be concluded 
that the high-z dust model can be discarded to a confidence level 
of 4<x for a homogeneous dust distribution. Also, it can be ob- 
served that the differences between the Q m o = 0.31, Qao = 0.69 
universe and the replenishing dust model are almost negligible. 
Finally, the differences between our model with respect to the 
Snla data are smaller than lcr. As discussed in §2, the advan- 
tage of the new model is that the dust distribution derives from 
current observational evidence. 

In Figure 4 the entire sample of photometric data from Snla 
is shown superposed on the theoretical models considered so 
far. It can be seen that the error bars of individual photomet- 
ric data include both the SFR-based dust and the 'replenishing' 
dust models. This model is quite flexible: the ambiguity in cur- 
rent cosmic star formation history studies (expulsion velocities, 
delay times and SFR curve parameterization) does not have a big 
impact on our results, especially if the constants are empirically 
fitted as explained in Sect. 2.3. 

In all these universe models (the grey dust ones), we have 
assumed a flat geometry with Q m o — 1 .0, but this assumption 
contradicts the results of 2dF (Peacock et al. 2001). For this 
reason we have also developed an open universe (£2 m o = 0.3, 
Qao = 0.0) with SFR-based dust whose differences are less than 
lcr compared to the Snla data (Figure 5 ). The implications of 
this open universe will be discussed in Sect. 4. 



Fig. 3. Averaged Snla data and theoretical models showing the modu- 
lus difference of distances with respect to a universe in constant expan- 
sion (Q. = 0). The high-z dust model is rejected at 4cr (dashed line). The 
best fit is given by an f2 mo = 0.31, Q. Ao = 0.69 universe (solid thin line) 
with no grey dust contribution. The SFR-based dust model represented 
takes into account a delay of 2.7 Gy between the dust formation in the 
disk and its homogeneous distribution in the IGM. 



3.3. Further evidence against high-z grey dust 

It might be argued that a non-homogeneous grey dust distri- 
bution could produce the deviation from the experimental data 
compared to Aguirre's model. According to this hypothesis, at 
high redshifts only the Snla whose line of sight to the observer 
passes through more voids (in other words: whose "beam" is 
emptier) can be observed, since the others are too absorbed by 
the dust. In this way, distant Snla are brighter than expected ac- 
cording to the high-z dust model, lying somewhere between the 
high-z grey dust and the £2 m o = 1 model, where the data points 
are located (Figure 4). 

To evaluate the feasibility of this hypothesis, the Snla sample 
is searched for close pairs of Snla with very different redshifts, 
in which the low z (z - 0.5) member of the pair is fainter than 
expected from the Milne model, while the high-z one (z > 0.9) 
is not. Were such a pair to be found, it should be concluded that 
inhomogeneities in the high-z grey dust model cannot be deemed 
responsible for the observed Snla distance modulus distribution. 

Assuming a Robertson-Walker metric, the angle subtended 
by an object in a flat universe is 



Dt(z) 
Dp(z)' 



(6) 



where Dj is the transversal distance and Dp the proper distance. 
From this equation, and assuming a cosmological model, the an- 
gle corresponding to a transversal distance can be evaluated. 

To minimize the influence of photometric and model errors 
when applying this test, it is convenient to look for apparent Snla 
pairs of very different redshifts, of which the Snla at low z show 
the highest brightness dimming. This maximum Snla dimming 
is observed around z = 0.5 (Riess et al. 1998, 2004; Perlmutter 
et al. 1999). Assuming a characteristic size of 30 Mpc for a void, 
the angle corresponding to this distance at z = 0.5 will represent 
the highest projected angular distance to search for a suitable 
pair to apply this test. This is a rather stringent requirement since 
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Fig. 4. Top: As in Figure 3, but representing the whole Snla sample. 
The SFR-based distribution fits the data within the error bars of indi- 
vidual data. Bottom: Detail of the < z < 0.7 region. 

a void size of 30 Mpc is smaller than those found by Einasto 
et al. (1994). From their data, an average diameter of around 
•frvoid - 91 Mpc is deduced. In Table 1, the angles corresponding 
to this distance at z = 0.5 for different models of universe are 
shown. 



£2m0 


£2ao 
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0.3 


0.0 


0.9502 


0.3 


0.7 


0.8451 


1.0 


0.0 


1.0155 



Table 1. Values for the angle (in degrees) subtended by a D T of 30 Mpc 
at z = 0.5 in different universes, with the assumed Hubble constant of 
65 km s -1 Mpc -1 . 



With these hypotheses, the requirements to apply this test to 
supernovae pairs are the following: 




0.0 0.5 1.0 1.5 2.0 

Z 



Fig. 5. Differences with respect to an empty universe. The flat (thick 
line) and open (dashed line) SFR-based dust universe models are within 
the lcr confidence level. 



1. The two members of the candidate pair have to be in the 
ranges 0.35 < z < 0.65 for the low z and z > 0.9 for the high 
z, respectively. This requirement maximizes the difference 
in brightness of the pair members to reduce the contribution 
from photometric errors as much as possible. 

2. For the same reason, the difference in redshift between the 
members of the pair must be Az > 0.4. 

3. The angular distance on the sky between the two members of 
the pair must be lower than the value of ^ 1 degree assumed 
in Table 1 . 

4. Both supernovae must be classified with the maximum con- 
fidence as type la. This requirement restricts the search to the 
Gold sample. 

For such a small Snla sample (157 Snla ) it might seem dif- 
ficult to fulfil these requirements. Luckily, the systematic search 
for high-z supernovae has centred on specific regions of the sky, 
mainly on the GOODS project HDFN and CDFS fields. The sky 
distribution of the Snla sample is represented in Figure 6. The 
whole sample is shown in the left panel, while on the right the 
Snla found in the HDFN are shown, demonstrating the small 
distances that separate the Snla in these deeply studied regions. 

Table 2 summarises the results of our study for different cos- 
mologies. Here, all the Snla at 0.35 < z < 0.65 that are members 
of a pair are weaker than expected in an empty universe. A rela- 
tively large number of pairs is found, considering that the whole 
sample is composed by only 186 Snla. 

4. Conclusions 

With the available Snla data we have calculated the best-fit dust- 
free universe model, obtaining Q m o = 0.31, Oao = 0.69. These 
results are very similar to those found in previous studies of lu- 
minosity distance to type la supernovae, leading to a reconcilia- 
tion of the age of globular clusters with the age of the universe 
(fo ^ 14.6 x 10 9 years, as derived from the model, when assum- 
ing a cosmological constant-like dark energy, cj = -1). Phantom 
energies can thus be avoided, as long as the mean values do not 
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Fig. 6. Location on the sky of the Snla sample. Left: Position of all the supernovae in the sample. Right: Detail of the Hubble Deep Field North 
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Table 2. Number of Snla pairs found with an angular distance between 
members lower than the angle 9 given by Table 1 for the corresponding 
cosmology. A: Pairs in which both members belong to the Gold sample. 
B: At least one of the members belong to the Silver sample. 



increase. However, the uncertainties in the age of globular clus- 
ters are still too great to constrain the equation of state of dark 
energy. 

A new grey dust model based on replenishing dust but in- 
cluding the observed comoving SFR density distribution has 
been derived, thereby providing an alternative model based on 
observational evidence. The extinction produced by this new 
SFR-based grey dust with Q mU = 1.0 is very close to the 
^mO = 0.31, £2ao = 0.69 model with no extinction. Apart from 
the difficulties in providing the mass required for a flat uni- 
verse and the age conflict arising therefrom, the models would be 
distinguished only by future supernovae surveys at high z with 
higher photometric accuracy, such as those to be carried up with 
SNAP {Supernova Acceleration Probe) (Aldering et al. 2002) 
and with the JWST (James Webb Space Telescope, Stockman 
et al. 1998), which could provide the evidence that will help to 
decide between a currently dark energy-dominated universe and 
a matter-dominated one. Moreover, the next generation Snla sur- 
veys, with more accurate data, will be essential because the un- 
certainties in the current Snla data, due to the contributions of 
approximations, photometric errors, luminosity maximum esti- 
mations at high z, etc., impose several restrictions on the validity 
of the cosmological results. Only when combined with other ex- 
periments such as CMB anisotropies or gravitational lenses is it 
possible to derive more stringent conclusions. 

All the flat grey dust models presented here consider 
Euclidean geometry with an adimensional mass parameter, 
£l m o — 1, a hypothesis against the data obtained so far (Peacock 
et al. 2001). Even with the results of del Popolo (2003), who 
has suggested that current Q m o values obtained from clusters of 



galaxies are underestimated by 20%, a parameter £2 m o — 1 is 
excluded from the acceptable margins. For this reason we have 
included an open model Q m o = 0.3, Qao = 0.0 with our SFR- 
based grey dust attenuation. Although the confidence level of 
this model with respect to the Snla data is below l<x, it intro- 
duces a serious conflict with the interpretation of the WMAP 
results (Spergel et al. 2003). To reconcile this model with the 
CMB anisotropies experiment it is necessary to resort to exotic 
mechanisms, such as thermalization of the radiation due to the 
intracumular dust at high z- 

The high-z dust model is rejected out to 4<x even if the 
grey dust presents a patchy distribution, as shown by the pan- 
test devised here. However, most of the data at very high red- 
shift corresponds to GOODS fields, which cover a small area 
on the sky (Figure 5). Pairs at other positions on the sky com- 
bined with high-accuracy photometry would allow us to ascer- 
tain whether there is any significant amount of grey dust in the 
universe. Nevertheless, the consistency of the standard model 
and the weak basis for the alternative explanations makes the 
possibility of having an open model describing our universe dif- 
ficult. 
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